Spin relaxation of two-dimensional electrons with a hierarchy of spin-orbit couplings 
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The density matrix formalism is applied to calculate the spin-relaxation time for two-dimensional 
systems with a hierarchy of spin-orbit couplings, such as Rashba-type, Dresselhaus-type and so on. 
It is found that the spin-relaxation time can be infinite if those coupling strengths a, (3, 71 and 72 
satisfy either condition (i) a = j3, 71 = or (ii) a = — j3, 72 = 0, which correspond to the vanishing 
Yang-Mills magnetic field. The effect caused by the application of an external magnetic field is also 
discussed. It is found that the longitudinal and in-plane spin components can possess infinite life 
time when the spin components, the Larmor precession frequency and the external magnetic field 
satisfy certain relations. 

PACS numbers: 72.25.Rb, 72.25.Dc 



I. INTRODUCTION 

Spintronics [lj, or spin based electronics |2|, has ab- 
sorbed much more attention during the last decade. One 
important issue in this research area is the manipulation 
of spin polarized electrons with the help of an electric 
field HI1I1I1I3 The system with spin-orbit couplings 
makes these efforts possible and thus brings great inter- 
ests from both academic and practical aspects 0- El S ■ 
However, the problem of the loss of the average micro- 
scopic spin is crucial in experimental data analysis and 
applicable device construction. The study of the spin- 
relaxation mechanisms of two-dimensional electrons is 
thus very important. 

Spin relaxation exhibits some properties of the spin 
dynamics, which plays inevitable role in realizing the ap- 
plicable spintronics devices. The main mechanism of spin 
relaxation in systems lacking inversion symmetry is the 
D'yakonov-Perel's mechanism |8|, |9j , in which the spin 
of the electron precesses due to an effective k-dependent 
magnetic field. For electrons in two-dimensional semicon- 
ductor heterostructures or quantum wells, the structure 
inversion as y mmetr y brings about the Rashba spin-orbit 
coupling Eamill, while the bulk inversion asymme- 
try in the A3B5 compounds leads to the Dresselhaus 01 
spin-orbit coupling. Spin-relaxation time in some semi- 
conductors with both Rashba and Dresselhaus couplings 
was calculated by analyzing the condition of spin de- 
cay [T3 . IT3 . and the effect of external magnetic fields was 
discussed |16| furthermore. An infinite spin-relaxation 
time was predicted in the system with equal Rashba 
and Dresselhaus coupling constants by making use of an 
SU(2) symmetry in fc-space. It is important to under- 
stand the spin-relaxation mechanism and the condition 
for infinite spin-relaxation time to occur, which would be 
helpful for overcoming the difficulties in the spin-based 
information processes. 

In this paper, we develop the aforementioned theory 
of spin relaxation to describe two-dimensional electron 
systems in the presence of U(l) Maxwell field and SU(2) 
Yang-Mills fields. Such a system can be realized in cer- 



tain semiconductor materials where the spin-orbit cou- 
plings, such as Rashba-type, Dresselhaus-type and etc., 
play crucial roles. Using density matrix formalism, we 
calculate the spin-relaxation time for the system with a 
hierarchy of spin-orbit couplings. In the absence of the 
Maxwell magnetic field, the infinity of the spin-relaxation 
time occurs if the spin-orbit couplings a, (3, 71 and 72 sat- 
isfy the condition in which the Yang-Mills magnetic field 
vanishes. In order to capture the physical essence of the 
emergence of an infinite spin-relaxation time, we further 
study the effect of the external magnetic field on the same 
systems and find that the longitudinal and in-plane spin 
components can also possess infinite life times when the 
spin orientation, the Larmor precession frequency and ex- 
ternal magnetic field satisfy some relations. Base on the 
analysis of spin-orbit systems with or without Maxwell 
magnetic fields, we expose a physics picture for a clear 
understanding of the infinite spin-relaxation time, which 
is helpful for the design of the spin-based devices. 

II. SPIN RELAXATION ARISING FROM 
SPIN-ORBIT COUPLINGS 

To start with a general formalism, we consider the 
Schrodinger equation for a particle moving in an external 
17(1) Maxwell field and an SU(2) Yang-Mills field 0, 

d 

ih—^(r,t) = H^(r,t), 

H = ^-(p - - r,A a f a ) 2 + eA Q + V A a r\ (1) 

where $ is a two-component wavefunction, A p = 
(A , Ai) denotes the vector potential of the Maxwell elec- 
tromagnetic field, and A M = A^T a that of Yang-Mills 
field with r a being the generators of SU(2) Lie group. 
It has been shown [Isj that the Yang-Mills fields can be 
realized in certain semiconductor materials. 

Firstly, we consider a two-dimensional system (in x- 
y plane) with a four-parameter Yang-Mills potentials: 

A = (0,0,0), A x = 2W0,/3 + a, 72 ), A y = 2f(/3 - 
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a,0, 71), ^ = (0, 0, 0) where a, (3 and 7 character- 
ize the strengths of spin-orbit couplings of Rashba-type, 
Dresselhaus-type, etc., respectively. If writing out the 
Hamiltonian explicitly, we have 



2 J.2 



H 



h 2 k 



+ V + k y <T x (a - j3) - k x a y (a + /3) 



2m 



21.2 



2m 



V 



with V = 2m(j 2 + 2(a< 
2(k y (a-(3),-k x [a + P), 



(2) 

7 2 2 )A 2 + v sc u n k = 

-72^))/^ and (J being 



the pauli matrix. Here m stands for the effective mass of 
electrons in the material, V sc = eA the scattering poten- 
tial which is independent of spin indices. The scattering 
is supposed to be elastic. The last term in above equa- 
tion H' = -|<7 • Jlk causes the precession of electron spins 
with the Larmor frequency fit which can be regarded as 
an effective magnetic field. 

To calculate the spin-relaxation time, we use density 
matrix formalism. The electron density matrix p(k) with 
components p ss '(k), s, s' be ing the indices of electron 
spin states, is denned by [T^.IT^| 



P(k) 



i[ff'(k), p(k)] + WWW - P(k')) = 0, (3) 



where r is the lifetime, Wkk' is the scattering probability 
from k to k' and the square bracket denotes commutator. 

Since H' contributes merely a small perturbation, 
the spin-relaxation time is much longer than the time 
for electron-momentum distribution to become isotropic, 
i.e., t Ti, n the momentum relaxation time. There- 
fore, it is convenient to split the density matrix into two 
parts, 



p = p + p' , with p' = 0. 



(4) 



Here we use a bar to denote an average taken over all 
directions of k and a prime to denote the deviation ones, 
p'(k) <C p- Taking average for Eq. ©, we have the fol- 
lowing relation, 



£ + ^[if'(k),p'(k)]=0. 



(5) 



Eq. can a l so be written out as 
P'(k) 



, -[ff'(k),p'(k)] - -[ff'(k),p'(k)] + ~[H'(k),p] 



(G) 



in which Eq. JSJ) has been used. Without taking account 
of the higher order terms, we need to solve the following 
equation 

UH'(k),p] + J2 WW//(k) - p'(k')] = 0. (7) 



This approximation is valid when Qki~i <C 1. For elastic 
process the scattering probability is a function of deflec- 
tion angle only, which makes it possible to expand the 
above equation in terms of Fourier series. After some al- 
gebra, one can express p' in terms of p. Substituting it 
into Eq. (J5J and employing the Boltzmann equation with 
only collision term, we can obtain the rate of averaged 
density matrix p, namely, 

1 -^r n [HL n , [H' n ,p]}, (8) 



dt ' sp - rel - 



h 2 



with 



K = f^H'(k)eM-^), (9) 
— = I dOWkv (1- cos n0), (10) 

T n J 

where cf>k is the angle between k and x-axis, and 8 — 

- 4>k> ■ 

Now we are in the position to investigate the kinetics 
of the spin density Si(t) = J <Zj(e, t)ds, where a,i(e,t) — 
[F + (e) — F_(e)]si(t), Si(t) = tT(a{p) and F±(e) refers to 
the distribution function projected along the direction 
parallel or anti-parallel to s = (sj, S2, S3). Accordingly, 
we obtain the evolution equation for the spin density at 
the time longer than t\ 15]: 

Si{t) = ~Sj{t) 

1 1 - J ds(F+ - F.)r n tx{[R'_ n , K,p]]°i} 



m 2h 2 ^ 



de{F+ -F-) 



(11) 

where i, j — x, y, z. 

For the Hamiltonian under consideration (0), we ob- 
tain the following, 

1 ^ 7i 2 +7 2 2 + ( a + /?) 2 

JL = 7i 2 +7 2 + (a-/?) 2 

T yy 
1 



-A. 
-A. 



= {a 2 +(3 2 )A 

7~ZZ 

1 1 



7i(a- 0) 



A. 



J_ _ J_ _ 72(Q + /3) 



A. 



(12) 



• xy 1 yx 

where the coefficient A is given by 



A = 



8m 



j de[F+(e) - F_(E)]n( E )e 



ft* 



de[F + (e) - F_(e)} 
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The above entities in Eq. I|12|) define the spin-relaxation 
tensor T = mat(^-). Diagonalizing this matrix (see Ap- 
pendix El for details), we obtain 



A 
2 



\ 



(13) 



with 
1 



Ti 



± 



ll,± 

- = 7 1 2 + 7l + 2(a 2 + /3 2 



\l (7i 2 + 7l) 2 + 8«/3( 72 2 - 7? + 2aP) 



(14) 



Clearly, two of the diagonal elements are always posi- 
tively definite and the other one 1/tii _ is not. The con- 
dition for a vanishing 1/tii turns out to be 

7l 2 (a + /3) 2 + 7 |(a-/3) 2 + (a + (3) 2 {a - (if = 0. (15) 
The above equation gives rise to two solutions 

(i) a = P, 7i = 0, 

(ii) a = -f3, 72 = 0. (16) 

Under these conditions, T~} and T~} are zero when the 

y y zz 

infinite spin-relaxation times emerge. 

Actually, the Yang-Mills "magnetic" field B z = b a r a 
can be calculated, namely 

4m 2 



6 2 = 



r) 2 h 2 
4m 2 

4m 2 



rph 2 



{P-a)j 2 , 
(a 2 -/? 2 ). 



(17) 



The condition i|15|) is also equivalent to 
|6| 2 = 0, b= (6i,6 2 ,6 3 ), 

which implies that the module of the Yang-Mills "mag- 
netic" field vanishes. In other words, the spin-relaxation 
time can be infinite when the module of the Yang-Mills 
magnetic field is null. This result is expected to be a 
criterion to evaluate whether there has an infinite spin- 
relaxation time in two-dimensional systems with spin- 
orbit couplings. In order to determine which spin com- 
ponent has an infinite life time, the spin precession needs 
to be analyzed concretely. 

The nonvanishing spin orbit coupling 71 or 72 will 
bring about some new features which may be useful for 
possible design with an infinite spin relaxation time. For 
the first case a = (3 and 71 = in Eq. I|ltj|) , the Hamilto- 
nian W reduces to 



H' 



-a ■ tth 



-k x (2aa y +720-2 



(18) 



The orientation of Larmor precession frequency Jlk is 
parallel to y'-axis as illustrated in Fig. |31 in appendix A, 
thus one can understand why the life time of the «Sy 
component is infinite while the other two decay. Here 
the y'-axis is defined in the diagonalization procedure of 
the spin-relaxation time tensor Ea. (|12f) for the first case 
given in appendix lAl 

From the appendix [S] and the figure therein, we can 
see that the strengths of spin-orbit coupling 72 and a 
determine the angle 9' between i~2k and z-axis. The angle 
9' can be manipulated by these two parameters, thus 
a definite alignment of spin with infinite life time can 
be realized with the help of tuning spin-orbit coupling 
strengths. On the other hand, the ratio of different type 
of spin-orbit coupling constants can be determined by 
means of measuring spin-relaxation time experimentally. 

For the second case in Eq. (|16|l . similar analysis can 
be carried out, which is omitted here. In the special case 
when both 71 and 72 vanish, one component of the tensor 
of spin-relaxation time becomes null, i.e., 1/t xx = or 
1/t to = 0, and thus the S x or S y has infinite life time, 
which is just the case considered in Ref. Q. 



III. SPIN RELAXATION AFFECTED BY 
EXTERNAL MAGNETIC FIELD 

In previous section, we considered the case of Aq = 
which means external magnetic field is absent. In the 
presence of the magnetic field, we should take account of 



A 



2VB 



(B X , By, B z ) 



(19) 



where fj,B is the Bohr magneton. The existence of an 
external magnetic field is known to affect the dynamics 
of the electron's spin. The Larmor precession of elec- 
tron's spin around a sufficiently strong longitudinal mag- 
netic field will suppress the precession about the internal 
random magnetic fields 20]. The cyclotron motion will 
change the wave vector k and affect the spin relaxation 
due to the DP mechanism. The density-matrix formal- 
ism is applicable for calculating electrons' spin-relaxation 
time. It will be convenient to expand the density matrix 
for electrons in terms of the unit and Pauli matrices, 



p(k) = / k + s k • a 



(20) 



where /k = ir[p(k)/2] is the spin-averaged-electron- 
distribution function and s k = ir[p(k)(?] is the spin per 
k-state electron. The kinetic equation for the spin distri- 
bution is given bv [11 l2ll E^] : 



dt 



s k x i^L + n k) + w c • [k x V k s k ] 

^2 Wkk'(Sk - Sk') = 0, 



(21) 



in which the second term refers to spin precession caused 
by spin-orbit couplings given in Eq. (J2J) together with 
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the external magnetic field; the third term is related to 
the wave vector variations due to the cyclotron motion 
whose frequency is u>c — and the last term denotes 
the collision integral. Since the internal random magnetic 
field is regarded as a perturbation, i.e., r2 k Ti <C 1, we can 
split the spin distribution function s k in the following, 



s° fc 



(22) 



where is a quasi-equilibrium distribution function and 
thus is independent of direction of k. Whereas, (5s k is 
a nonequilibrium correction arising from spin-orbit cou- 
plings as well as other internal random magnetic fields 
and thus it contains only first angular harmonics of the 
spin distribution [23j because elastic scattering processes 
is taken into account only, accordingly 



(5s k = Ri cos(4> k ) + R-2 sin(0 fc ), 



(23) 



where the two vectors Ri and R2 are irrelevant to the 
direction of the wave vector k though they are functions 
of the module of k in general. Substituting Eq. 1)221) and 
(|23)l into Eq. Ij21[l. we obtain the following equations: 



dt 
d5s k 

dt 



+ s2 x u) L + <Ss k x fi k = 0, 



s° x f2 k + 5s k x u) L 



uo c ■ [k x V k <5s k ] 



fe k 

n 



0. 



(24) 



(25) 



where t\ is momentum-relaxation time whose definition 
is also given by Eq. I|1U|) for n = 1. In the light of the 
number of total electrons N = 2 ^ k fa and the single 



electron spin S° = 



E 



N 



summing Eq. I)24|) over the 



wave vectors, we obtain the balance equation describing 
electron's spin relaxation 



dS° 
~dT 



S° x uj l + f S° = 0, 



(26) 



where the spin-relaxation tensor f referring to the inverse 
of spin-relaxation times is defined as 



(27) 



The nonequilibrium correction 5s k can be obtained from 
Eq. (|25|l in which the contribution of the rate dSs^/dt 
is negligible because its magnitude is of higher order in 

f2 k Tl. 

Firstly, we rotate the original coordinate {£, y, z} 
which is related to the principal crystal axes to the new 
one {x',y',z'} (illustrated in Fig. The coordinates in 
both systems are related, (x,y,z) = (£', y' , z')R T , by 



cost/ cos ip, — smiy9, sinfcos^ 
R = I cos 9 sin (p, cos ip, sin 9 sin <p 
— sin#, 0, cos# 



(28) 



Here, 9 is the angle between z and z', and <p the angle be- 
tween y' and y([010]). Similar equations are valid for mo- 
mentum components, ki — Rijk'j (here i,j = x,y,z). It 
is convenient to calculate the nonequilibrium correction 
<5s k and the components of the spin-relaxation tensor in 
the new frame of coordinate where the Larmor frequency 



vector u>l 



L Bjij in the original coordinate becomes 



ujl = ujlz' in the new coordinate. 




FIG. 1: The scheme of coordinate frames, z-axis parallel 
with the [001] growth axis. 6 and cp are the polar and the 
azimuthal angles of the external magnetic field B. z 1 is chosen 
in alignment with the orientation of B, y' is lying in x-y plane, 
and x 1 is chosen to form a right-hand triple with y 1 and z' . 

After tedious calculation, we obtain the spin-relaxation 
tensor (inverse of the spin-relaxation time) T for degen- 
erate electrons with Fermi energy Ep which is given in 
appendix^] These results are valid for arbitrary random 
internal magnetic field and arbitrary orientation of the 
external field, from which we obtain several conclusions 
that will be illustrated respectively. 



A. Longitudinal relaxation 

The longitudinal spin-relaxation rate is 1/t l — T z i z i. 
From Eq. HB1(I . we can obtain the following conclusion: 
T z i z i — when either 

a = j3, 71 = 0, tp = tt/2, 9 = tan" 1 (2a ! /72), 




Schematic of the case 



0. 



P, 7i = 0, 



Fig. 

tp — tt/2, 8 — tan 1 (— 72/20?), the external magnetic field B 
is parallel to ^'-direction, and x' axis is antiparallel to the 
Larmor frequency fl k . 
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or 

a = -{3, 72 = 0, if = 0, 0= -tsaT 1 {2a hi). (29) 

When the longitudinal spin-relaxation time t l is infinite, 
the spin component S z > has an infinite life time. Cer- 
tainly, the Hamiltonian H' describing the electron spin 
precession arising from spin-orbit coupling can also be 
written as Eq. I)18|l for the former case in Eq. 1)29(1 . One 
can see from Fig. © that S z > is the component parallel 
to fit when 6' = 9 = tan -1 (20/72) and cp = it/ 2. Thus 
the infinite life time of S z > can easily be understood from 
physical point of view, that is to say, S z > will not pre- 
cess about Larmor precession frequency ilk and external 
magnetic B when Jlk is parallel with B. While the other 
components S x i and S y > have finite life time due to the 
precession around l~ik arising from internal random mag- 
netic field. If both of the strength of spin-orbit couplings 
71 and 72 vanish, the longitudinal component of the ten- 
sor T z i z i is zero when 9 = tt/2, tp = tt/2, which can be 
seen from the Eq. 1B1|) . This recovers the special case 
discussed in Ref. jlq . 

For the latter case in Eq. (|29|l . the Hamiltonian H' 
becomes 

H' = -a ■ f2k = k y (2aa x - 71 o>). (30) 

2k 

The Larmor frequency fik = ^^(2a, 0, —71) 

h 

parallels to the external magnetic field B((9 = 
— tan _1 (2a/7i), ip = 0). So S z > does not decay 
because Jlk parallels to B. 

B. Transverse relaxation 

Let us analyze the spin relaxation in the plane perpen- 
dicular to the external magnetic field. One can find that 
the transverse components of spin-relaxation tensor can 
also be zero (i.e., T x > x > — 0, T x > y > — 0, Y y i x i — 0) when 
either 

a = /3, 7 i =0,(^ = V2,^ = -tan- 1 (^), 

or 

a=-/3 ! 7 2 = 0,^ = 0,e = tan- 1 (^). (31) 

la 

As shown in Fig. |3 the x' axis is antiparallel to the 
Larmor frequency (i.e., x' \\ — l~2k) for the former case 
in Eq. I|31fl . The spin component S x > does not precess 
about the Larmor precession frequency ilk- We know 
that B is perpendicular to Ilk from Fig. |3 thus S x > will 
precess about the constant external magnetic field B in 
the plane paralleling with ilk- So the random internal 
magnetic fields and external magnetic field can not in- 
duce spin relaxation for the spin component S x i. Then 
S x i has an infinite life time associating with T x i x i = 0. 
And the admixture of the x' component to y' component 



of spin-relaxation tensor (which is described by Y x i y i , 
T y i x i) are zero, namely T x i y i = and Y y i x i = which 
can also be calculated from Eqs. (|B1IB5|) . 

For the latter case in Eq. I|31|) , the direction of the ex- 
ternal magnetic field is also perpendicular to the Lar- 
mor frequency ilk as illustrated in Eq. (|3U|I . The spin 
component S x > is antiparallel to the Larmor frequency 
(x 1 || — ilk). The random internal magnetic field and ex- 
ternal magnetic field will not induce spin relaxation for 
the spin component S x > due to the same reason as afore- 
mentioned. Then S x > has an infinite life time, and the 
components of spin-relaxation tensor T x t y i, T y ' x i vanish 
(see Eqs. (IB1IB50 . 

The components T x > z > and T y i z i are smaller than oth- 
ers when the external magnetic field is sufficiently strong 
(f^Ti <C wl). Under this condition the in-plane spin 
components rapidly rotate and the admixture of the in 
plane components to z'-component (which is described by 
r x 'z' , Tj/'z') plays no role in the spin dynamics. There- 
fore the above result manifests the general solutions of 
spin-relaxation time for the Hamiltonian Eq. J3J| . 



IV. SUMMARY 

In the above, we developed a consistent theory of spin 
dynamics to describe particles moving in an external 
[7(1) Maxwell field and an SU{2) Yang-Mills field which 
characterizes spin-orbit couplings in certain semiconduc- 
tors (such as Rashba-type, Dresselhaus-type or other 
complex types). We used the density-matrix formalism 
to calculate the spin-relaxation time in such systems in 
the absence and in the presence of the external magnetic 
field, respectively. In the absence of external magnetic 
field, we find that the spin component S y > or S z i has 
an infinite life time if the strengths of spin orbit cou- 
plings a, P, 71 and 72 satisfy either (i) a = f3, 71 = 
or (ii) a — —f3, 72 = 0. In such a case, the Yang-Mills 
magnetic field vanishes. From these conditions, the di- 
rection of the spin component with infinite life time can 
be manipulated by tuning a, 72 or /3, 71 respectively. In 
the presence of the external magnetic field, we consid- 
ered the magnetic effect on the two-dimensional system. 
We obtained that the longitudinal spin-relaxation time 
is infinite when the S Z ' is parallel to ilk and B if cither 
(i') a = /3, 71 = 0, ip = tt/2, 9 = tan- 1 (2a/7 2 ) or (ii') 
a = —(3, 72 = 0, (p = 0, 9 = — tan _1 (2a/7i). By making 
use of analysis in detail, we conclude that the in-plane 
spin component S x i can also have infinite life time if S x ' 
is antiparallel with f2k and perpendicular to B for either 
(iii') a = 0, 71 = 0, ip = tt/2, 6 = tan" 1 (-7 2 /2a) or (iv') 
ol = -P, 72 = 0, cp = 0, 6 = tan" 1 (71 /2a). That is to 
say, electron spins have infinite life times if they precess 
in the plane parallel to the vector of the Larmor preces- 
sion frequency arising from spin-orbit couplings. These 
solutions provide a better understanding on the spin dy- 
namics of two-dimensional system with a four-parameter 
Yang-Mills potentials, which characterizes a hierarchy of 
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spin-orbit coupling in certain semiconductor materials. 
It is expected to expose some more clues for manipulat- 
ing spin via certain spin-orbit couplings in semiconduc- 
tors or elaborating spintronics storage devices with long 
spin-relaxation time. 
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Hence, the new spin components S — U S is obtained 




APPENDIX A: ON THE DIAGONALIZING 
BASES 



We can easily obtain the matrix U 
/ 1 





u = 



2a 



72 







V^f + 4a 2 yj-fl + 4a 2 
72 2a 



(Al) 



V V 72 2 + 4a 2 y/j* + 4a 2 / 



which diagonalizes the spin-relaxation tensor given by 
Eq. H2) to be Eq. JT3J|, 



T = U^T U 



(A2) 



This matrix U turns the evolution equation for the spin 
density Ea. lfTTJl to be, 



dS' _ v -idS_ 
dt ~ dt 

= — T _1 S" 



-(U- 1 TU)U- 1 S 



(A3) 



Fig. 3: The orientation of the Larmor frequency f2k for the 
case a = /3, 71 = 0, which defines the y'-axis of the new frame 
of coordinate 



S' y = smd'S y + co$e'S z , 

S' z = -cos9'S y + sm6'S z , 
where tan 9' = 2a This means the existence of spin- 
orbit coupling 72 makes the orientation of the Larmor 
frequency fik to change from y to y'-axis. As illustrated 
in Fig. |21 the 9' refers to the angle between the Larmor 
frequency and z-axis. 



APPENDIX B: THE SPIN-RELAXATION 
TENSOR 



In the presence of magnetic held, we obtained the fol- 
lowing spin-relaxation tensor 



2fc 2 n 
h?D+D- 



-{[1 + t 2 (uj 2 + u>l)] x [{a + (3) 2 cos 2 <p + (a - (3) 2 sin 2 ip 



+ (72 sing — (a + (3) cos9 simp) 2 + (7! sing + (a — j3) cos 9 cos ip) 2 ] 
+4tIujc<jJl[{oi 2 - P 2 ) cos 9 + (a + /3)7i singcosy + {fi - a)^ 2 sin 9 sin ip]} 



(Bl) 



2fc 2 n 1 + {u; 2 + lu 2 l )t 2 



r, 2 



-{ 



x [(72 cos 9 + (a + 0) sin 9 sin ip) 2 



, n , n \ ■ n n2i (a — (3) 2 sin 2 ip + (a -\- [3) 2 cos 2 (p 
+ (71 cosg- (a- f3) sing cosy) 2 ] + ^ — — 2 V 2 ; -} 



(B2) 



r„ v = + ^ x [(72 cos + (a + /?)singsiny) 2 + ( 7l cosg - (a - (3) sing cosy) 2 



1 



1 + LU 2 C T 2 



x [(72 sinff — (a + (3) cos sin ip) 2 + (71 sing + (a — (3) cos cosy) 2 ]} 



(B3) 



T x , y , = ^P~{— ~jj~J} ^ ULTl x [(72 cos 9 + (a + P) sin 9 sin ip) 2 + (71 cos fl - (a - 0) sin fl cos y) 2 ] 

H ^5 — 5- x [—(a 2 — (3 2 ) cosOuicTi + ((a — j3) 2 — (a + (3) 2 ) cos 6 cos u> sin u> 

1 + ui z c t{ 

+71 [(a — /3) sin # sin <^ — (a + [3)luc t i sintfcos^] + 72 [(a — 0)^c T \ sin sin ^ + (a + (3) sin 9 cos ip]]} 

(B4) 

ry x , = [1 - UJ ^ T ^ UJLTl x [( 72 cos 9 + (a + (3) sin 6> sin cp) 2 + ( 7l cos fl - (a - /3) sin 9 cos y) 2 ] 

x [(a 2 - f3 2 ) cos ^wcri + ((a - /3) 2 - (a + f3) 2 ) cos cos <p sin 



1 + uj 2 c tI 



+7i[(a — /3) sin^sin^ + [a + $)wct\ sin 9 cos </j] + 72 [(a — /3)wc' r i sin sin ip + (a + 0) sin cos if]]} 

(B5) 



r 

with fc = \J2ti\Ef and £>± = 1 + (ul ± c^c) 2 ''" 2 ■ 
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